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Abstract 

We study symplectic deformations of Gabor frames, using the co- 
variance properties of the Heisenberg operators. This allows us to 
recover in a very simple way known results. We thereafter propose a 
general deformation scheme by Hamiltonian isotopies, which are paths 
of Hamiltonian flows. We define and study in detail a weak notion of 
Hamiltonian deformations, using ideas from semiclassical analysis due 
to Heller and Hagedorn. This method can be easily implemented using 
symplectic integrators. 

1 Introduction 

The theory of Gabor frames (or Weyl-Heisenberg frames as they are also 
cahed) is a rich and expanding topic of harmonic analysis. It has many 
applications in time- frequency analysis, signal theory, and mathematical 
physics. The aim of this article is to initiate a systematic study of the 
symplectic transformation properties of Gabor frames, both in the linear 
and nonlinear cases. Strangely enough, the use of symplectic techniques 
in the theory of Gabor frames is very often ignored; one example (among 
many others) being Casazza's seminal paper [9] on modern tools for Weyl- 
Heisenberg frame theory, where the word "symplectic" does not appear a 
single time in the 127 pages of this paper! There are however exceptions: 
in Grochenig's treatise [27j the metaplectic representation is used to study 
various symmetries; the same applies to the recent paper by Pfander et al. 
[45], elaborating on earlier work [29] by Han and Wang, where symplectic 
transformations are exploited to study various properties of Gabor frames. 

In this paper we consider the notion of deformation of a Gabor system 
using the Hamiltonian isotopies we introduce in section [5l A Hamiltonian 
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isotopy is a curve {ft)o<t<i of difFeomorphisms of phase space M^"" starting 
at the identity, and such that there exists a Hamiltonian function H, usuahy 
time-dependent such that the (generahzed) phase flow {ft^)t determined by 
the Hamilton equations 

X = dpH{x,p,t) , p = -dpH{x,p,t) (1) 

consists of precisely the mappings ft for < t < 1. It follows that a Hamilto- 
nian isotopy consists of symplectomorphisms (or canonical transformations, 
as they are called in Physics). Given a Gabor system ^((/>, A) with window 
(or atom) (j) and lattice A we want to find a working definition of the defor- 
mation of ^(0, A) by a Hamiltonian isotopy {ft)o<t<i- While it is clear that 
the deformed lattice should be the image A^ = /t(A) of the original lattice 
A, it is less clear what the deformation <j)t = fti(p) of the window <j) should 
be. A clue is given by the linear case: assume that the mappings ft are 
linear, i.e. symplectic matrices St; assume in addition that there exists an 
infinitesimal symplectic transformation X such that St = e^^ for < t < 1. 
Then {St)t is the flow determined by the Hamiltonian function 

H{x,p) = -^{x,pfjX{x,p) (2) 

where J is the standard symplectic matrix. There exists a one-parameter 
group of unitary operators {St)t satisfying the operator Schrodinger equation 

= H{x, -ihdx)St 

where the formally self-adjoint operator H[x, —ihdx) is obtained by replac- 
ing formally p with —ihdx in ©; the matrices St and the operators St 
correspond to each other via the metaplectic representation of the symplec- 
tic group. This suggests that we define the deformation of the initial window 
(j) hy (j)t = St<j3. It turns out that this definition is satisfactory, because it 
allows to recover, setting t = 1, known results on the image of Gabor frames 
by linear symplectic transformations. This example is thus a good guideline; 
however one encounters difficulties as soon as one want to extend it to more 
general situations. While it is "reasonably" easy to see what one should do 
when the Hamiltonian isotopy consists of an arbitrary path of symplectic 
matrices (this will be done in section U]), it is not clear at all what a "good" 
definition should be in the general nonlinear case: this is discussed in section 
[6l where we suggest that a natural choice would be to extend the linear case 
by requiring that (j)t should be the solution of the Schrodinger equation 



associated with the Hamiltonian function H determined by the equahty 
{ft)o<t<i = {fi^)o<t<i', the Hamiltonian operator H would then be associ- 
ated with the function H by using, for instance, the Weyl correspondence. 
Since the method seems to be difficult to study theoretically and to imple- 
ment numerically, we propose what we call a notion of weak deformation, 
where the exact definition of the transformation (j) i — > (pt of the window 
(j) is replaced with a correspondence used in semiclassical mechanics, and 
which consists in propagating the "center" of a sufficiently sharply peaked 
initial window (j) (for instance a coherent state, or a more general Gaussian) 
along the Hamiltonian trajectory. This definition coincides with the defini- 
tion already given in the linear case, and has the advantage of being easily 
computable using the method of symplectic integrators (which we review 
in section 15. 3p since all what is needed is the knowledge of the phase flow 
determined by a certain Hamiltonian function. Finally we discuss possible 
extensions of our method. 

Notation and terminology 

The generic point of the phase space M^" = M" x M" is denoted hy z = {x,p) 
where we have set x = (xi, p = {pi, ...,pn). The scalar product of 

two vectors, say p and x, is denoted by p • x or simply px. When matrix 
calculations are performed, z,x,p are viewed as column vectors. We will 
equip M?"^ with the standard symplectic structure 



the n X n zero and identity matrices). The symplectic group of M^" is 
denoted by Sp(n); it consists of all linear automorphisms of M^" such that 
a{Sz,Sz') = a{z,z') for all z,z' £ M?^. Working in the canonical basis 
Sp(n) is identified with the group of all real 2n x 2n matrices S such that 
S'^JS = J (or, equivalently, SJS^ = J). 

We will write dz = dxdp where dx = dxi • • • dx^ and dp = dpi ■ ■ ■ dpn- 
The scalar product on L^(M") is defined by 



and the associated norm is denoted by 1 1 • 1 1 . The Schwartz space of rapidly 
decreasing functions is denoted by 5(M") and its dual (the space of tempered 
distributions) by 5'(M"). 



a[z, z ) = p ■ X — p • x; 



in matrix notation a{z,z') = {z')'^Jz where 
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2 Gab or Frames 



Gabor frames are a generalization of the usual notion of basis; see for in- 
stance Grochenig [S^, Feichtinger and Grochenig [16], Balan et al. [3], Heil 
[32], Casazza [9] for a detailed treatment of this topic. In what follows we 
give a slightly modified version of the usual definition, better adapted to the 
study of symplectic symmetries. 

2.1 Definition 

Let 4>he a non-zero square integrable function (hereafter called window) on 
M", and a lattice A in M^n^ 2.e. a discrete subset of M^". The associated 
?i-Gabor system is the set of square-integrable functions 

g(0,A) = {f''(z)0:zG A} 

where T^{z) = ^-^'^i^^^)/^ jg the Heisenberg operator. The action of this 
operator is explicitly given by the formula 

(see e.g. [221 IMl [39] ; it will be justified in section g^])- We will cah the 
Gabor system Q{g, A) a /i-frame for L^(]R"), if there exist constants a, 6 > 
(the frame hounds) such that 

a||^|P< j;i(V|f'^(zo)0)|'<fe||^|P (4) 

20 e A 

for every square integrable function tjj on M". When a = b the /i- frame 
^(^f, A) is said to be tight. 

Remark 1 The product {'ijj\T^{zQ)<j)) is, up to the factor {27rh)~"', Wood- 
ward's cross-ambiguity function I50j : it is up to a (symplectic Fourier trans- 
form) the cross-Wigner distribution W{'tp,<p) as was already observed by 
Klauder see IM\M^- 

2.2 Rescaling properties 

For the choice h = 1/2tt the notion of ?i-Gabor frame coincides with the 
usual notion of Gabor frame as found in the literature. In fact, in this case, 
writing T{z) = T^l'^'^ [z] and p = uj, have 

mf{z)<i^)\ = mT{z)<p)\ 
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where t{z) is the modulation operator defined by 

for zo = {xo,ujq). The two following elementary results can be used to go 
from one definition to the other: 

Proposition 2 Let = i^r. n '^^r )• '^he system Q{(I),K) is a Gabor 

yU ZTTfil J 

frame if and only if Q{(f),D^A) is a h- Gabor frame. 

Proof. We have f''{xo,27rhpo) = f (xo,Po) where f (xo,Po) = f^/'^''{xo,Po)- 
By definition Q{(f>,A) is a Gabor frame if and only if 

20 e A 



for every tp G L^(M") that is 



(a;o,po)eA 



this inequality is equivalent to 



a||V||'< E mTHxo,2Trhpo)(l>)f<bM\' 

{a:o,Po)eA 



that is to 



aiivii'< E Mf^{xo,pom'<b\m' 

(a;o,(27rS)-ipo)GA 

hence the result since {xq, {2Trh)~^po) G A means that {xo,po) G D^A. m 

We can also rescale simultaneously the lattice and the window ( "change 
of Planck's constant"): 

Proposition 3 Let Q((p,A) be a Gabor system, and set 

/(x) = {2Trh)-''/^^{x/V2^). (5) 
Then Q{(f),A) is a frame if and only if Q{(l)^, s/^TrhA) is a h-frame. 
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Proof. We have = -^1/^/2;;^/ o*^ where -^1/^/2^/ ^ Mp(n) has projec- 
tion 

on Sp(n) (see Appendix A). The Gabor system Q{(j)^, \/27r/iA) is a ?i-frame 
if and only 

aiivii' < E Mn^^)A/v2^ni,M^ < 

zoeVS-Tr/iA 

for every ^ £ L^(M"'), that is, taking the symplectic covariance formula ([6|) 
into account, if and only if 

a\m' < Yl l(M^,_o^|f((2vr;i)-V2^o,(2vr;.)V2p,)^)|2 < 5||^||2. 
But this is inequality is equivalent to 

and one concludes using Proposition O ■ 

Remark 4 In Appendix A we state a rescaling property for the covering 
projection vr'^ : S 1 — > S of metaplectic group Mp(n) onto Sp(n) (formula 

m)- 

3 Symplectic Covariance 

The following formula will play a fundamental role in our study of symplec- 
tic covariance properties of frames. It relates Heisenberg-Weyl operators, 
linear symplectic transformations, and metaplectic operators (we refer to 
Appendix A for a concise review of the metaplectic group Mp(n) and its 
properties). 

Let S G Mp(n) have projection tt^'XS) = S € Sp(n). Then 

f^\z)S = Sf^\S^'^z). (6) 

For a proof see e.g. \22\ [Ml [39] : one easy way is to prove this formula 
separately for each generator J, M^^m, Vp of the metaplectic group. 
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3.1 A first covariance result 



Gabor frames behave well under symplectic transformations of the lattice 
(or, equivalently, under metaplectic transformations of the window). Let 
S £ Mp(n) have projection S G Sp(n). The following result is well-known, 
and appears in many places in the literature (see e.g. Grochenig [27], Pfander 
et al. [35]). Our proof is somewhat simpler since it exploits the symplectic 
covariance property of the Heisenberg-Weyl operators. 

Proposition 5 Let (/) £ L'^{R"-) (or (p G 5(IR");. A Gabor system g{<p,A) 
is a h- frame if and only if Q{S(p,SK) is a h- frame; when this is the case 
both frames have the same hounds. In particular, Q{(j),K) is a tight h- frame 
if and only ifQ{(j),A) is. 

Proof. Using the symplectic covariance formula ([6]) we have 

z&SA zeSA 

and hence, since Q{(p,A) is a ?i- frame, 

a||§-VlP< E mTHz)S^)\^<b\\S-'ij\\^ 
zeSA 

The result follows since ||S'~"^VII = IIV'II because metaplectic operators are 
unitary; the case (p £ 5(1^"") is similar since metaplectic operators are linear 
automorphisms of 5 (M"). ■ 

Remark 6 The result above still holds when one assumes that the window (p 
belongs to the Feichtinger algebra So{W^) (see Appendix B and the discussion 
at the end of the paper). 

3.2 Gaussian frames 

The problem of constructing Gabor frames Q{(p,K) in I/^(]R) with an ar- 
bitrary window (j) and lattice A is difficult and has been tackled by many 
authors (see for instance the comments in |28j . also |45j). Very little is 
known about the existence of frames in the general case. We however have 
the following characterization of Gaussian frames which extends a classical 
result of Lyubarskii [H] and Seip and Wallsten ^\ : 
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Proposition 7 Let (j)Q{x) = (7r/i)~"/^e~l^l Z^'' (the standard centered Gaus- 
sian) and Aai3 = oX^ x /SZ" with a = (qi, a„) cLnd (3 = (/3i, T/ien 
^(00, Aa/s) is a frame if and only if ajPj < 2T:h for 1 < j < n. 

Proof. Bourouihiya [5j proves this for h = l/27r; the result for arbitrary 
h> follows using Proposition [3l ■ 

It turns out that using the result above one can construct infinitely many 
symplectic Gaussian frames using the theory of metaplectic operators: 

Proposition 8 Let(pQ be the standard Gaussian. The Gabor system Q{(()q, K(^p) 
is a frame if and only if Q{S(j)Q, SA^fs) is a frame (with same bounds) for 

every S G Mp(n). Writing S in block-matrix form ( ^ ^ ) the window 



S(f)Q is the Gaussian 



Sct>^^^x) = (l^Y'^ (detX)i/4e-2k(^+^>')-- (7) 

where 

X = -{CA^ + DB^){AA^ + BB^)-^ (8) 

Y = [AA'^ + BB^)-^ (9) 

are symmetric matrices, and X > 

Proof. That a((/'o,Aa/3) is a frame if and only if ^(S'(/>q, S'A^^) is a frame 
follows from Proposition [5l To calculate S<j)Q it suffices to apply formulas 
([6i]) and ([65]) in Appendix A. ■ 



3.3 An example 

Let us choose h = 1 /2tt and consider the rotations 

S^^f-ost sin A 
y— smr cost J 

(we assume n = 1). The matrices (St) form a one-parameter subgroup of 
the symplectic group Sp(l). To {St) corresponds a unique one-parameter 
subgroup (St) of the metaplectic group Mp(l) such that St = 7r^/^'^(S't). It 
follows from formula ()Aip in Appendix A that Stcj) is given for t ^ kn (k 
integer) by the explicit formula 

J — CO 
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where m{t) is an integer (the "Maslov index") and 



W(x, x',t) = (ix^ + x'^) cos t - 2xx') . 

2sint ^ 

Remark 9 The metaplectic operators St are the "fractional Fourier trans- 
forms" familiar from time-frequency analysis (see e.g. Almeida /i/, Namias 

Applying Proposition [5] we recover without any calculation the results of 
Kaiser [35] (Theorem 1 and Corollary 2) about rotations of Gabor frames; 
in our notation: 

Corollary 10 Let G{(j), A) be a frame; then Q{S(p, SA) is a frame for every 
5 G Mp(n). 

4 Symplectic Deformations of Gabor Frames 

The symplectic covariance property of Gabor frames studied above can be 
interpreted as a first result on Hamiltonian deformations of frames because, 
as we will see, every symplectic matrix is the value of the flow (at some 
time t) of a Hamiltonian function which is a homogeneous quadratic poly- 
nomial (with time-depending coefficients) in the coordinates Xj,pk. We will 
in fact extend this result to deformations by afhne flows corresponding to 
the case where the Hamiltonian is an arbitrary quadratic function of these 
coordinates. 

4.1 The linear case 

The first example in subsection l3.3l (the fractional Fourier transform) can be 
interpreted as a statement about continuous deformations of Gabor frames. 
For instance, assume that St = e*^, X in the Lie algebra sp(n) of the 
symplectic group Sp(n) (it is the algebra of all 2n x 2n matrices X such 
that XJ + JX^ = 0; when n = 1 it reduces to the condition Tr X = 0; see 
[19\ I22j). It is then easy to check that 5^ can be identified with the flow 
determined by the Hamilton equations z = JdzH for the function 

H{z) = -\z^{JX)z (11) 

that is 

^St = XSt. (12) 
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A fundamental observation is now that to the path of symplectic matrices 

t I — 7- St, < t < 1 corresponds a unique path t i — St, < t < 1, of 
metaplectic operators such that Sq = I and Si = S (see Appendix A). This 
path satisfies the operator Schrodinger equation 

ih^^St = HSt (13) 

where H is the Weyl quantization of the function H. Collecting these 
facts, one sees that Q{S(j)Q, SA^is) is obtained from the initial Gabor frame 
Q{^Q,Acii3) by a smooth deformation 

t^g{St4,StAafi) , o<t<i. (14) 

More generally, let S be an arbitrary element of the symplectic group Sp(n). 
The latter is connected so there exists a path (in fact infinitely many) 
t I — > St, < t < 1, joining the identity to S in Sp(n). An essential result, 
generalizing the observations above, is the following: 

Proposition 11 Let t i — > St, < t < 1, be a path in Sp(n) such that 
Sq = I and Si = S. There exists a Hamiltonian function H = H{z,t) such 
that St is the phase flow determined by the Hamilton equations z = JdzH. 
Writing 

the Hamiltonian function is the quadratic form 

H{z,t) = \{btCj -CtDj)x^ 

+ {DtAj - CtBj)p ■x + \ {BtAj - AtBj)p' (16) 

where At = dAt/dt, etc. 

Proof. The matrices St being symplectic we have Sj JSt = J ■ Differ- 
entiating both side of this equality with respect to t and setting we get 
SfJSt + SfJSt = or, equivalently, 

jstsr' = -{sj^r'sj^j = (jStsrY. 

This equality can be rewritten JStSf^ = (JStS^^)^ hence the matrix 
JStS^^ is symmetric. Set JStS^'^ = Mt = Mf; then 

St = XtSt , Xt = -JMt (17) 
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(it reduces to formula ()12p when Mt is time-independent). Define now 

H{z,t) = -\z^{JXt)z- 

using (jl7p one verifies that the phase flow determined by H consists precisely 
of the symplectic matrices St and that H is given by formula (|16|) . ■ 

Remark 12 Formula \10^) also follows from the more general formula i31\) 
about Hamiltonian isotopies in Proposition U5\) below. 

Exactly as above, to the path of symplectic matrices t i — > St corresponds 
a path t I — > St of metaplectic operators such that Sq = I and Si = S 
satisfying the Schrodinger equation (I13p . Thus, it makes sense to consider 
smooth deformations (jl4p for arbitrary symplectic paths. This situation will 
be generalized to the nonlinear case in a moment. 

4.2 Translations of Gabor systems 

A particular simple example of transformation is that of the translations 
T{zo) : z I — > z -\- zq in R^". On the operator level they correspond to 
the Heisenberg-Weyl operators T^{zq). This correspondence is very easy to 
understand in terms of "quantization" : for fixed zq consider the Hamiltonian 
function 

H{z) = a{z, Zq) =p-xo-po- X. 

The corresponding Hamilton equations are just x = xq, p = po whose solu- 
tions are x(t) = x(0) + txo and p{t) = p{0) + tpo, that is z{t) = T{tzo)z{0). 
Let now 

H = a{z, Zq) = {-ihdx) -xo-po-x. 
be the "quantization" of H, and consider the Schrodinger equation 

ihdt(j) = a-(z, zo)(l). 

Its solution is given by 

(l){x,t) = e-*'"(^'^°)/V(^, 0) = f''\tzo)(l)ix,0) 

(the second equality can be verified by a direct calculation, or using the 
Campbell-Hausdorff formula [H [22l IMl [39]). 

Translations act in a particularly simple way on Gabor frames; we write 
T{zi)A = A + zi. 
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Proposition 13 Let zq,zi G M^". A Gabor system Q{cl),A) is a h-frame if 
and only if g{T^{zo)(p,T{zi)A) is a h-frame; the frame bounds are in this 
case the same for all values of zo,zi. 

Proof. We will need the following well-known [191 1121 [Ml ES] properties of 
the Heisenberg-Weyl operators: 

f^\z)f^{z') = e'''^''^'^/^'f^{z')f^'{z) (18) 
f^z + z') = e-'''^'^'''^''^^f^\z)f^{z'). (19) 

Assume first zi = and let us prove that Q{T^''{zQ)(f), A) is a ?i-frame if and 
only Q{(j),K) is. We have, using formula ([T8|) and the unitarity of T^\zq), 

Y,mT\z)f\z,m^ = |(V|e^'^(^'^")/'^f '"X^o)T''(z),/.)| 

= Y,mT\zo)f\zm 

= Y,\iTH-zoWHzm; 

it follows that the inequality 

a\m' <Y;^mf\z)f^{zom' <bm\^ 

zeA 

is equivalent to 

aM\'<Y^Mf\z)^)\'<bUf 

hence our claim in the case zi = 0. We next assume that zq = 0; we have, 
using this time formula (jl9p . 

zeT{zi)A z6A 

= YmTHz,)f\zM' 

z6A 

= Y,\iTH-zi)^\f'\zm' 

zGA 

and one concludes as in the case zi =0. The case of arbitrary zq,zi imme- 
diately follows. ■ 



12 



Identifying the group of translations with M^" the inhomogeneous (or 
affine) symplectic group ISp(n) [8l[19] is the semi-direct product Sp(n) xM^"; 
the group law is given by 

(S, z){S' , z') = {SS' , z + Sz'). 

Using the conjugation relation 

S-'T{zo)S = T{S-'zo) (20) 

one checks that ISp(n) is isomorphic to the group of all affine transformations 
of of the type ST{zo) (or T{zo)S) where S G Sp(n). 

The group ISp(n) appears in a natural way when one considers Hamil- 
tonians of the type 

H{z,t) = \M{t)z- z + m{t)- z (21) 

where M{t) is symmetric and m{t) is a vector. In fact, the phase flow 
determined by the Hamilton equations for (j2ip consists of elements of ISp(n). 
Assume for instance that the coefficients M and m are time-independent; 
the solution of Hamilton's equations i = JMz + Jm is 

z^ = e'J^^zo + (JM)-i(e*'^*^ - I)JM (22) 

provided that det M / 0. When det M = the solution (p2|) is still formally 
valid and depends on the nilpotency degree of X = JM. Since X = JM G 
5p{n) we have St = e*^ G Sp(n); setting = X~^(e*^ - I)u the flow (//^) 
is thus given by 

fi" = T{^t)St G ISp(n). 

The metaplectic group Mp(n) is a unitary representation of the double 
cover Sp2(n) of Sp(n) (see Appendix A). There is an analogue when Sp(n) 
is replaced with ISp(n): it is the Weyl-metaplectic group WMp(n), which 
consists of all products T{zo)S; notice that formula ([6]), which we can rewrite 

S-^f^'{z)S = f^iS-h). (23) 

is the operator version of formula (I20p above. 

5 The Group Ham(n) 

In this section we review the basics of the modern theory of Hamiltonian 
mechanics from the symplectic point of view; for details we refer to [12\ [331 
H6] ; we have also given elementary accounts in [221 . 
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5.1 Hamiltonian flows: properties 

A Hamiltonian system ([T]) can be written in compact form as 

z = Jd^H{z,t) (24) 

where J is the standard symplectic matrix. The Hamiltonian function H is 
assumed to be twice continuously differentiable in z, and continuous in t. 
We denote by the mapping R^" — > R^" which to an initial condition 
zo associates the value z = //^(-Zo) of the solution to (plj) at time t. The 
family {ft^)t of all these mappings is called the phase flow determined by 
the Hamiltonian system (I24j) . 

It is often useful to replace the notion of flow as defined as above by that 
of time- dependent flow (fHi)'- fHi is the function such that flli{z') is the 
solution of Hamilton's equations with z{t') = z' . Obviously 

/5' = /fo(A^o)"' = /f U^)"' (25) 
and the f^, satisfy the groupoid property 

fi!t'fe,t' = fi!t" ' fi!t = (26) 

for all t, t' and t" . Notice that it follows in particular that {fl\>)~^ = f^f 

An essential property which links Hamiltonian dynamics to symplectic 
geometry is that each mapping //^ is a diffeomorphism such that 

[Dfi'{z)fjDfi'{z) = Dfl'{z)J [Df^{z)f = J. (27) 

Here Dfl^{z) is the Jacobian matrix of the diffeomorphism //^ calculated 
at the point z = {x,p): 

dz ~ d{x,p) 

if zt = fl^{z). The equality (I27p means that the matrix DfP{z) is sym- 
plectic: DfP{z) G Sp(n) for every z and t. Any diffeomorphism / of phase 
space R^" satisfying the condition 

Df{zfjDf{z) = Df{z)JDf{zf = J (28) 



is called a symplectomorphism. Formula (|27p thus says that Hamiltonian 
flows consist of symplectomorphisms, which is a well-known property from 
classical mechanics [21 [Ml US] • 

A remarkable fact is that composition and inversion of Hamiltonian flows 
also yield Hamiltonian flows: 



14 



Proposition 14 Let {ft) and {ft) he the phase flows determined by two 
Hamiltonian functions H = H{z,t) and K = K{z,t). We have 

fl'ff = ff*'' w^th H#K{z,t) = H{z,t) + K{{ft'')-\z),t). (29) 
(/f)~' = /f ^^th H{z,t) = -H{ff{z),t). (30) 

Proof. It is based on the transformation properties of tlie Hamiltonian fields 
Xh = JdzH under diffeomorphisms; see [221 ES SS] for detailed proofs. ■ 
We notice that even if H and K are time-independent Hamiltonians, 
then H^K and H are time-dependent. 



5.2 Hamiltonian isotopies 

Formula (I28p above shows, using the chain rule, that the symplectomor- 
phisms of M^" form a group, which we will denote by Symp(n). 

Let us now focus on the Hamiltonian case. We will call a symplecto- 
morphism / such that / = //^ for some Hamiltonian function H and time 
t = 1 a Hamiltonian symplectomorphism. The choice of time t = 1 in this 
definition is of course arbitrary, and can be replaced with any other choice 
t = a: we have / = fl^" where Ha{z, t) = aH{z, at). 

Hamiltonian symplectomorphisms form a subgroup Ham(n) of the group 
Symp(n) of all symplectomorphisms; it is in fact a normal subgroup of 
Symp(n) as follows from the formula 

g-'fl'a = /f °^ 

valid for every symplectomorphism 5 of M^" [Ml ESI El]. That Ham(n) 
really is a group follows from the two formulas ()29p and pop in Proposition 
[T4l above. 

The following result is, in spite of its simplicity, a deep statement about 
the structure of the group Ham(n). It says that every continuous path of 
Hamiltonian transformations passing through the identity is itself the phase 
flow determined by a certain Hamiltonian function. 

Proposition 15 Let {ft)t be a smooth one-parameter family of Hamiltonian 
transformations such that /o = Id- There exists a Hamiltonian function 
H = H{z,t) such that ft = //^. More precisely, {ft)t is the phase flow 
determined by the Hamiltonian function 

H{z, t) = -J^ z'^J {ft o ft-') (Az)dA (31) 

where ft =dft/dt. 
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We refer to Wang [l9] who gives an elementary proof of formula (|3ip. 
The result goes back to the seminal paper of Banyaga [2], but the idea is 
already present in Arnold [2] (p. 269) who uses the apparatus of generating 
functions. 

We will a call smooth path {ft)t in Ham(n) joining the identity to some 
element / G Ham(n) a Hamiltonian isotopy. 

5.3 Symplectic algorithms 

Symplectic integrators are designed for the numerical solution of Hamilton's 
equations; they are algorithms which preserve the symplectic character of 
Hamiltonian flows. The literature on the topic is immense; a well-cited 
paper is Channel and Scovel [52j. Among many recent contributions, a 
highlight is the recent treatise [36j by Kang Feng, Mengzhao Qin; also see 
the comprehensive paper by Xue-Shen Liu et al. |51j, and Marsden's online 
lecture notes [42] (Chapter 9). 

Let (//^) be a Hamiltonian flow; we assume first that H is time-independent 
so that we have the one-parameter group property fl^f^ = f^tt- Choose 
an initial value zq at time t = 0. A mapping /a* on M^" is an algorithm 
with time step-size At for (//^) if we have 

fKti^) = fAt{z) + 0{At''); 

the number k (usually an integer > 1) is called the order of the algorithm. In 
the theory of Hamiltonian systems one requires that /a* be a symplectomor- 
phism; /^t is then called a symplectic integrator. One of the basic properties 
one is interested in is convergence: setting At = t/N {N an integer) when 
do we have lim.N^ooift/N)'^ i^) = ft^i^)'^ One important requirement is 
stability, i.e. {ft/N)^ i^) must remain close to z for small t (see Chorin et 
al. [10]). 

Here are two elementary examples of symplectic integrators. We assume 
that the Hamiltonian H has the physical form 

Hix,p) = Uip) + Vix). 
• First order algorithm. One defines {xk+i,Pk+i) = fAt{xk,Pk) by 

Xk+i = Xk + dpU{pk - dxV{xk)At)At 
Pk+i =Pk- d^V{xk)At. 
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• Second order algorithm. Setting 

x'k = Xk + \dpU{pk) 

we take 

Xk+i = Xf. + \dpU{pk) 
Pk+i = Pk - dxV{x',,)At. 

One can show, using Proposition [15] (Wang [l9]), that both schemes are 
not only symplectic, but also Hamiltonian. For instance, for the first order 
algorithm above, we have f^t = /a* where K is the now time-dependent 
Hamiltonian 

K{x,p,t) = U{p) + V{x-dpU{p)t). (32) 

When the Hamiltonian H is itself time-dependent its flow does no longer 
enjoy the group property fi^f^ = f^t'^ redefine the notion of 

algorithm in some way. This can be done by considering the time-dependent 
flow (/^^,) defined by ([25]): /^^, = //^ One then uses the following 

trick: define the suspended flow (//^) by the formula 

f?{z',t') = {fli,,,{z'),t + t'); (33) 

one verifies that the mappings (the "extended 

phase space") satisfy the one-parameter group law fl^f^ = f^^, and one 

may then define a notion of algorithm approximating //^ (see Struckmeier 
[l8j for the extended phase space approach). 

6 Hamiltonian Deformations of Gabor Systems 

Let / G Ham(n) and (/t)o<t<i a Hamiltonian isotopy joining the identity 
to /; in view of Proposition [15] there exists a Hamiltonian function H such 
that ft = ft^ for < t < 1. We want to study the deformation of a /i-Gabor 
frame g{(j),A) by (/t)o<t<i; that is we want to define a deformation 

gi<p,A)^giUt^,ftAy, (34) 

here Ut is an (unknown) operator associated in some (yet unknown) way 
with ft- The fact that when ft = St £ Sp(n) we have 

g{(t>,A) ^g{St4>,StAy, (35) 

where St G Mp(n) with St = ir'^iSt) suggests that: 
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• The operators Ut should be unitary; 

• The deformation should reduce to 135\) when the isotopy {ft)o<t<i 
lies in Sp(n). 

The following property of the metaplectic representation gives us a clue. 
Let (St) be a Hamiltonian isotopy in Sp(n) C Ham(n). We have seen in 
Proposition [TT] that there exists a Hamiltonian function 

H{z,t) = \M{t)z-z 

with associated phase flow precisely {St)- Consider now the Schrodinger 
equation 

ih^ = H^l, , ^(•,0)=V'o 

where H is the Weyl quantization of H (it is a formally self-adjoint operator 
since H is real). It is well-known [22l [231 IS] that = •S'tV'o where {St)t is 
the unique path in Mp(n) passing through the identity and covering {St). 
This suggests that we should choose {Ut)t in the following way: let H be the 
Hamiltonian function determined by the Hamiltonian isotopy {ft)'- ft = ft^ ■ 
Then quantize H into a operator H using the Weyl correspondence, and let 
{Ut) be the solution of Schrodinger's equation 

in^Ut = HUt , Fo = /d. 
at 

The operators Ut are unitary. Let in fact u{t) = {Util^\Utip) where ■0 is in 
the domain of H (we assume it contains 5(M")); we have 

ihu{t) = {HUtMUtij) - {UtMHUtij) = 

since H is (formally) self-adjoint; it follows that {Utip\Utip) = {ip\ip)- 

While definition ()35p of a Hamiltonian deformation of a Gabor system 
is "reasonable", its practical implementation is difficult because it requires 
the solution of a Schrodinger equation. We will try to find a weaker, more 
tractable definition of the correspondence (3^, which is easier to implement 
numerically. 

6.1 The semiclassical approach 

The idea of this method comes from semiclassical mechanics; historically 
it seems to be due to Heller ^33j . Hagedorn j30l EI] has given the method 
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a firm mathematical basis; also see Littlejohn's review paper [39j and the 
construction in |44j . 

For fixed zq we set zt = [zq) and define the new Hamilton function 

H,,{z,t) = {d,H){zut){z-zt) + \DlH{zut){z-ztf- (36) 

it is the Taylor series of H at zt with terms of order and > 2 suppressed. 
The corresponding Hamilton equations are 

i = Jd,H{zt,t) + JDlH{zt,t){z - Zt). (37) 

We make the following obvious but essential observation: in view of the 
uniqueness theorem for the solutions of Hamilton's equations, the solution 
of (|37p with initial value zq is the same as that of the Hamiltonian system 

i{t) = Jd,H{z{t),t) (38) 

with z(0) = Zq. Denoting by (/^ ^'') the Hamiltonian flow determined by 

we thus have //^(zq) = //^^"(zq)- More generally, the fiows {ff^°) and 
(//^) are related by a simple formula: 

Proposition 16 The solutions of Hamilton's equations [3l\ ) and (E^j ^'"s 
related by the formula 

zit) = Zt + StiziO) - zo) (39) 

where zt = f^ [zq), z{t) = fl^ [z) and {St)t is the phase flow determined by 
the quadratic time- dependent Hamiltonian 

H^{z, t) = \DlH{zt, t)z ■ z. (40) 

Equivalently, 

fl'{z) = T[zt-St{zQ)]St{zm (41) 
where T{-) is the translation operator. 

Proof. Let us set u = z — zt. We have, taking (I37p into account, 

n + it = Jd,H{z{t),t) + JDlH{zt,t)u 
that is, since zt = JdzH{zt,t), 

ii = JDlH{zt,t)u. 
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It follows that u{t) = St{u{0)) and hence 

z{t) = fl'izo) + Stu{0) = zt- Stizo) + St{z{0)) 

which is precisely (j39|) . ■ 

The nearby orbit method (at order = 0) consists in making the Ansatz 
that the approximate solution to Schrodinger's equation 

where 

< = f\zo)4 (42) 
is the standard coherent state centered at zq is given by the formula 

^(x,t) = e^(*'^«)f^(zi)^t(zo)T'"X^o)-V^o (43) 
where the phase 7(t, zq) is the symmetrized action 

l{t,zo)= f {\a{zt,,zt,)-H{zt,,t'))dt' (44) 

JO 

calculated along the Hamiltonian trajectory leading from zq at time to = 
to zt at time t. One shows that under suitable conditions on the Hamiltonian 
H the approximate solution satisfies, for |t| < T, an estimate of the type 

m;t)-i^{-M<C{z^,T)^\t\ (45) 

where C{zq,T) is a positive constant depending only on the initial point zq 
and the time interval [— T, T] (Hagedorn \2>Q\ I31j). 

Remark 17 Formula shows that the solution of Schrodinger's equation 
with initial datum (pQ is approximately the Gaussian obtained by propagating 
(/>Q along the Hamiltonian trajectory starting from z = while deforming it 
using the metaplectic lift of the linearized flow around this point. 

6.2 Application to Gabor frames 

Let us state and prove the main results of this paper. 

In what follows we consider a Gaussian Gabor system G{4>oj ^)j ^-Pplying 
the nearby orbit method to 0q yields the approximation 

^ ei^mf\zt)St<i)^^ (46) 

where we have set St = St{0). Let us consider the Gabor system ^((/)^,Aj) 
where Aj = //^(A). 
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Proposition 18 The Gabor system Q{(f)^,A.t) is a Gabor h-frame if and 
only ifQ{(j)Q,A) is a a Gabor h-frame; when this is the case both frames have 
the same bounds. 

Proof. Writing 

ItW = E mTHz)4)\' (47) 

z&At 

we set out to show that the inequahty 

2 ^ r /„i.\ ^ U\ L/.l |2 



<im <hm\' (48) 

(for all ■0 £ holds for every t if and only if it holds for t = (for all 

^ G L^(]R")). In view of definition (j46p we have 

/t(V') = Ei(^i^'(^)^'(^*)^*'^o)h 

the commutation formula (jlSp yields 

f^\z)f^{zt) = e''''-''''^/^f''\zt)f''\z) 

and hence 

= 5^|(V|f'^(^t)f''(/,^(z))5,0S^)|2. 
Since T^{zt) is unitary the inequality ()48p is thus equivalent to 

a\m' < mf\fi'{z))St4r < b\m'. m 

In view of formula ([5U|) we have, since StZQ = because zq = 0, 

fi'iz) = Stz + fi'{0) = Stz + zt 
hence the inequality ()49p can be written 

a||V'|P< J]|(V^|f'"X5t^ + ^t)5t</.S^)|'<6||V'll'. (50) 
zeA 

In view of the product formula (fT9]) for Heisenberg-Weyl operators we have 

f^iStz + zt) = e'''^^''^''^/'^''f^\zt)f''\Stz) 
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so that ([50|) becomes 

all^lP < Y.mf\z,)f\Stz)St4)\'' < ftllV'lh (51) 

the unitarity of T^^{zt) implies that ()5ip is equivalent to 

aUf < Y.mf\S,z)St<l^'^\'' < hM\\ (52) 

Using the symplectic covariance formula ^ we have 

f^{Stz)St = Stf\z) 
so that the inequality (|52p can be written 

aim' <Y,mStfHz)4)\' <bM\'; 

zGA 

since St is unitary, this is equivalent to 

The Proposition follows. ■ 

The fact that we assumed that the window is the centered coherent 
state (pQ is not essential. For instance, Proposition 1131 shows that the result 
remains valid if we replace (pQ with a coherent state having arbitrary center, 
for instance (j)^^^ = T^{zQ)(f)Q. More generally: 

Corollary 19 Let Q{(p,A) be a Gabor system where the window (p is the 
Gaussian 

where M = , ImM > 0. Then ^(0[S A*) is a Gabor h-frame if and only 
if it is the case for Q{(l), A). 

Proof. It follows from the properties of the action of the metaplectic group 
on Gaussians (see Appendix A) that there exists S G Mp(n) such that = 
S4>Q. Let 5* = vr''(5) be the projection on Sp(n) of S; the Gabor system 
g((/>|f,A) is a n-frame if and only if ^(5" Vl/, -^^^A) = g{(j)^^, S'^ A) is a 
h-hame in view of Proposition [5l The result now follows from Proposition 

m ■ 
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7 Discussion and Additional Remarks 



We have given one working definition of tlie notion of Hamiltonian defor- 
mation of a Gabor frame; this definition uses ideas from semiclassical me- 
chanics. However, we have used no approximations. We could therefore call 
this deformation scheme "weak Hamiltonian deformation". An important 
remark is that in all our results, one can assume that the window (p belongs 
to the Feichtinger algebra S'o(M"') (reviewed in Appendix B). This is due to 
the fact that we have transformed the Gabor frames under consideration 
only by the phase space translations T^^{z) and by metaplectic operators; it 
turns out that S'o(M"') is the smallest Banach algebra invariant under these 
operations, and thus semiclassical propagation preserve the Feichtinger al- 
gebra (see de Gosson [23J). A consequence is that the weak Hamiltonian 
deformation scheme behaves well with respect to the Feichtinger algebra. It 
is unknown whether this property is conserved under passage to the general 
definition ([35]) . that is 



where Ut is the solution of the Schrodinger equation associated with the 
Hamiltonian operator corresponding to the Hamiltonian isotopy {ft)o<t<i- 
This because one does not know at the time of writing if the solution to 
Schrodinger equations with initial data in 5o(IR") also is in 50(1^") for given 
time t. 

Since our definition of weak deformations was motivated by semiclassical 
considerations one could perhaps consider refinements of this method using 
the asymptotic expansions of Hagedorn [30l[3l] and his followers; this could 
then lead to "higher order" weak deformations, depending on the number 
of terms that are retained. Still, there remains the question of the general 
definition (fM|) where the exact quantum propagator is used. It would indeed 
be more intellectually (and also probably practically!) satisfying to study 
this definition in detail. As we said, we preferred in this first approach to 
consider a weaker version because it is relatively easy to implement numeri- 
cally using symplectic integrators. The general case (j54p is challenging, but 
not probably out of reach. From a theoretical point of view, it amounts to 
construct an extension of the metaplectic representation in the non-linear 
case; that such a representation indeed exists has been shown in our paper 
with Hiley j25j (a caveat: one sometimes finds in the physical literature a 
claim following which such an extension could not be constructed; a famous 
theorem of Groenewold and Van Hove being invoked. This is merely a mis- 
understanding of this theorem, which only claims that there is no way to 




(54) 



23 



extend the metaplectic representation so that the Dirac correspondence be- 
tween Poisson brackets and commutators is preserved). There remains the 
problem of how one could prove that the deformation scheme (j54p preserves 
the frame property; a possible approach could consisting in using a time- 
slicing (as one does for symplectic integrators); this would possibly also lead 
to some insight on whether the Feichtinger algebra is preserved by general 
quantum evolution. 

APPENDIX A: METAPLECTIC GROUP AND GAUS- 
SIANS 

Let Mp(n) be the metaplectic representation of the symplectic group Sp(n) 
(see [ini [22l [23] ) ; it is a unitary representation of the double cover Sp2 (^^) 
of Sp(n): we have a short exact sequence 

— ^ Z2 — > Mp(n) ^ Sp(n) 

where vr'* : S 1 — > S is the covering projection; we explain the appearance of 
the subscript H below. The metaplectic group is generated by the following 
elementary operators: 

• Fourier transform: 

J^={^T" I e--'/^^(x')dx' (55) 

(notice the presence of the imaginary unit i in the prefactor); 

• Unitary dilations: 

M^^„^ = i"»y/|detL|V'(La;) (detL^O) (56) 

where m is an integer depending on the sign of det L: m G {0, 2} if 
det L > and m G {1, 3} if det L < 0; 

• ''Chirps": 

Vp = e-'^'-^'/^^'tPix) (P = P^). (57) 

The projections on Sp(n) of these operators are given by 7r'^(J) = J, 
n'''iML,m) = Ml, and ■k''\Vp) = Vp with 

(the matrices Vp are sometimes called "symplectic shears"). 
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The projection of a covering group onto its base group is defined only up 
to conjugation; our choice -and notation- is here dictated by the fact that to 
the /i-dependent operators ([55]) and ([56]) should correspond the symplectic 
matrices (j58p . For instance, in time- frequency analysis it is customary to 
make the choice h = l/27r. Following formula relates the projections tt^ and 
vr = ^1/2-: 

TTiS) = AMyV2^nSM^) (59) 

wliere = 

Metaplectic operators are not only unitary operators on L^(M") but also 
linear automorphisms of which extend by duality to automorphisms 

of 5'(M"). 

There is an alternative way to describe the metaplectic group Mp(n). 

Let 

W{x, x') = \Px^ -Lx-x + iQx"^ (60) 

where P,L,Q are real n x n matrices, P and Q symmetric and L invertible 
(we are writing Px^ for Px ■ x, etc.). Let m be a choice of argdetL as in 
formula (j56p : each S € Mp(n) is the product to two operators of the type 

Sw,mHx) = {^^^^'^"'^Ad^\ [ e-*^(^'^')/'^V(^')dx' (Al) 

(see [Ml [221 [21]). The operators Sw, 

m can be factorized as 

Sw,m = V.pML,mJW.Q (A2) 

and hence belong to Mp(n). The projection Sw = '^^{Sw,m) is characterized 
by the condition 



(x,p) = Sw{x,p') 



p = dxW{x, x') 
p' = —dxW{x, x') 



this condition identifies W with the generating function of first type, familiar 
from Hamiltonian mechanics [21 [22l [Ml [38] . A straightforward calculation 
using the expression (I60p of W yields the symplectic matrix 

PL-^Q-L^ PL-\ 



Sw 



The metaplectic group acts on Gaussian functions in a particularly sim- 
ple way. Let M be a complex n x n matrix; we assume in fact that M 
belongs to the Siegel half-space 

S+ = {M : M = M'^,ImM > 0}. 
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We call generalized centered coherent state a Gaussian function of the type 
and for zq G M?^ we set 

A,.o = THzo)4j (62) 

(it is a Gaussian centered at the point zq). The symplectic group Sp(n) acts 
transitively on the Siegel half-space via the law |19) 

{S, M) ^ a{S)M = {C + DM){A + BM)-^. (63) 

One can show p3] that if M = X + then 

X = -{CA^ + DB^){AA^ + BB^)-^ (64) 
Y = {AA^ + BB^)-^. (65) 

This action induces in turn a transitive action 

of the metaplectic group Mp(n) on the set Gn of Gaussians of the type 
These actions make the following diagram 

Mp(n) X G„ Gn 

i i 
Sp(n)xS+ ^ S+ 

commutative (the vertical arrows being the mappings {S,4'\,i ^q) ' — ^ {Si^) 
and <?5>|/^p i — > M, respectively). 

The formulas above can be proven by using either the properties of the 
Wigner transform, or by a calculation of Gaussian integrals using the oper- 
ators Sw,m defined by formula ()Aip . 

APPENDIX B: FEICHTINGER'S ALGEBRA 

The Feichtinger algebra S'o(M"') was introduced in [13 1 114 1 [T5]: it is an impor- 
tant particular case of the modulation spaces defined by the same author; 
we refer to Grochenig's treatise [27j for a complete study of these impor- 
tant functional spaces. Also see Feichtinger and Luef [18j for an up to date 
concise review. 
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The Feichtinger algebra is usuahy defined in terms of short-time Fourier 
transform 

V^xPiz) = I e-2-'f^V(x')(/'(x' - x)dx'; (66) 
which is related to the cross- Wigner transform by the formula 




where ^^/27rl('^) ~ il^{xV2TTh) and (f)^ (x) = (f)(—x); equivalently 

VM^) = {^r'"^-'^'-''W{^y^,A\/^){z^). (68) 

The Feichtinger algebra 5o(M") consists of all ■0 S 5'(M") such that 
V^ip G L^(]R^") for every window (f). In view of the relations (|67p . (j68p 
this condition is equivalent to VF('i/',0) S L^iM?"-). A function G -L^(M") 
belongs to 5'o(M") if and only if W^) G ^^(M^"); here = W{il),^) is the 
usual Wigner function. The number 

\mls, = W{i^,'i>)\\L^M.^r.)= j \W{i^,^){z)\dz (69) 

is the norm of 0" relative to the window (j). We have the inclusions 

5(M") C 5o(M") c C7°(M") n ^^(IR") n F{L^{W)) (70) 

where F(Li(M")) is the image of L^(M") by the Fourier transform. One 
proves that 5*0 (M") is an algebra, both for pointwise multiplication and for 
convolution. 

An essential property of the Feichtinger algebra is that it is closed under 
the action of the Weyl-metaplectic group WMp(n): if S 5*0 (IR"), S G 
Mp(n), and zq e W we have both 50 G Sq{W) and T^^zq)^ G 50(1^"). In 
particular G Sq{W) if and only if G 50(1^"). One proves that Sq{W) 
is the smallest Banach space containing 5(M"') and having this property. 
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